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FINDING USABLE PORTION OF SIGMOID CURVE 

CROSS REFERENCE TO RELATED APPLICAllON 
This application claims the benefit of U.S. Provisional Patent Application No. 60/456,613 
5 filed March 21, 2003» which is incorporated herein by reference. 

TECHNICAL FIELD 
The technical field relates to analyzing observations of a sample to determine the 
concentration of a substance within the sanq>le» and more particularly to determining a concentration 
10 via a sigmoid curve representing observations for the sample. 

BACKGROUND 

One way of determining the concentration of a substance in a sanq)le is by performing serial 

dilution on the san^le. Serial dilution techniques collect a finite number of data points for the sample 
15 by taking one or more observations (e.g., indicating optical density) of various dilutions (e.g., 

dilutions formed by adding various quantity of dihient to the sample). For exanqple, dilutions of 10%, 

1%, 0. 1 %, etc. can be measured for optical density. 

The results can then be used to detennine a concentration of the substance in the sample via 

reference to a sigmoid curve representing serial-dilution observations for a sample having a known 
20 conc^itration of the substance (sometknes called a "standard" or "characteristic" sigmoid curve). 

FIG. 1 shows such a sigmoid curve 120. The sigmoid curve can be represented by the four-paramet^ 

Formula (1). 



The parameters of Formula (1) can be chosen so that the function / (x) calculates the optical density 
25 based on a particular dilution x. Given an optical density 130 for a sample having an unknown 
concentration of the substance and fhs degree of dilution associated with the sample, the 
concentration of the substance can be back-calculated. In practice, plural observations of the optical 
density can be taken for plural degrees of dilution and applied to ib& standard curve. 

30 Various techniques have been \ised to define the curve, analyze the observations, and 

calculate a concentration. One method is described by O'Connell, et al., "CahTsration and assay 
development using the four-parameter logistic model," Chetnometrics and Intelligent Laboratory 
Systems, 20 (1993) 97-1 14, Elsevier Science Publishers B-V., Amsterdam C'O'Conneir), which is 
hereby incorporated herein by reference. The O'Connell approach describes determining a rrrnihnmmi 

35 detectable conc«itration (MDC) and a reliable detection limit (RDL). 
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The O'CamicU technique can produce significant variability in its results. In certain 
scenarios, variability is to bt avoided. llierefore» there exists a need for technologies that avoid 
variability in their results and otherwise make better use of the sigmoid curve. 

5 sui^imahy 

Various technologies described herein proceed by analyzing data via a sigmoid curve, such 
as a standard sigmoid curve representing observations of a sample having a known quantity of a 
substance. 

Observation of a test sample can be taken (e.g., in serial dilution scenarios) for a sample, and 

0 the technologies described herein can be applied when calculating concentration of a substance 
within the sam})]e via the sigmoid curve» 

In some examples, derivative techniques are applied to detcnnine various points or ranges of 
the curve. The points or ranges can represent the usable portion of the sigmoid curve or a threshold 
value reliably indicating presence of the substance. 

5 Some of the technologies exhibit superior reduction in variability of results and are thus 

useful in a variety of fields, such as testing for levels of antibodies or antigens in a sample. For 
example, the technologies can be used when testing for titers of antibodies or antigens (e.g., in serum) 
via serial dilutioa The technologies can also be applied with advantage to a variety of scenarios, 
such as bioassay, polymerase chain reactive (TCR'O assay, radioimmuno assay ("RIA*')> cell grov^ 

10 assay, cell death assay, enzyme-linked immunosorbent assay (**£LISA'*), toxin neutralization assay 
("TNA") (e.g., scenarios involving detenmning the concentration of biological toxins such as 
anthrax), and flow cytometry scenarios, or any test in which san^le results are derived from a 
standard sigmoid curve. 

Additional features and advantages of the disclosed technologies will be noade apparent from 

15 the following detailed description of illustrated emibodimcnts, which proceeds with reference to flie 
accon^anying drawings. 

BBIEF DESCKIFnON OF THE DRAWINGS 

FIG. 1 shows an exemplary sigmoid curve. 
10 FIG. 2 is a flowchart of an exenylaiy method for deteTmining concenliation in a test saiiq)le 

via a sigmoid curve, 

FIG. 3 is a flowchart of an exeizq)lary method for determining the usable portion of a 
sigmoid curve. 

HG. 4 shows an exenq>laiy implen3entation of a method for determining the usable portion 
\5 of a sigmoid curve, such as tibie method shown in FIG. 3. 

FIG. 5 shows an cxsmplary implementatian of a method for determining a threshold titer on 
a sigmoid curve. 
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HG. 6 is a flowchart of an exemplary method for using derivatives to determine the 
concentration of a substance in a test sample. 

FIG. 7 shows an exenQ)lary sigmoid curve employed in a back calculation technique. 

FIG. 8 shows an exemplary sigmoid curve fit to data. 
5 FIG. 9 shows an exemplary curve showing titer computation. 

HG. 10 shows an exemplary plot illustrating back-cailculation of concentrations. 

HG. 11 is a block diagram showing an exeinplary system for carrying out the technologies 
described herein. 

FIG. 12 is a screen shot of an exen^lary user interface for presenting results of the 
1 0 technologies described herein. 

FIG. 13 is a screen shot of an exenqplaiy user interface for presenting results of tib^ 
technologies as applied to a standard sigmoid curve. 

FIG. 14 is a screen shot of an exemplary user inter&ce for presenting results of the 
technologies as applied to a test sigmoid curve. 
15 FIGS. 15andl6areascrecnshotofanexcnq)laryuserinterfeceforpresentingresuitsofthe 

technologies in text form. 

DETAELED DESCRIPTION 
Example 1 - Exemplary Overview 

20 In various of the exanq>le$ described herein, calculations for detennimng concentration of a 

substance in a test sample are performed by identifying a usable portion of a signooid curve. 
Performing the calculations with the usable portion of the curve can improve or optimize the 
concentration calculation. The portion of the curve that is usable can be determined via derivative. 
FIG. 2 shows an exemplary mediod 200 for determining the concentration of a substance in a test 

25 ssraple via a sigmoid curve. The method 200 can be in^lemented in software. 

At 210, the usable portion of a standard sigmoid curve is determined. For exan^le, a range 
of points between two endpoints can be defined as tiie usable portion. The standard sigmoid curve 
can be generated via observations of a sample having a known anaount of the concentcaticm present 
Observations can be taken for different dihitions (e.g., dilution ratios). In practice, phnal 

30 observations can be taken for each dilution. To fecilitate calculations related to the curve, a curve 
represented via four parameters {e.g., such as those of Formula 1) can be fit to the observations. 

At 220, the concentration (e.g., amount or quantity) of a substance in a test sample can be 
determined via fte usable portion of the curve. For exanqjle, given a set of observations of the test 
san^>le, some of the observations may relate to portions of the curve outside of the usable area. Such 

35 observations need not be included in concentration calculation. Instead, a subset of the observations 
(e.g., those associated with the usable portion of the curve) can be used in the calculations to 
determine concentration. 
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Example 2 - Exemplary Determination of Usable Portion 
The usable portion of a sigmoid curve (e.g., a standard sigmoid curve)* can be determined 
via derivatives of the sigmoid curve. For example, a maximum and minimum of the aecoud 
derivative can be used. FIG. 3 shows an exemplary method 300 by vvhich the usable portion of a 
5 sigmoid curve can be determined. The method 300 can be implemented in software. 

At 310, a first bound (e.g, endpoint) of the range is found via the second derivative of tlic 
sigmoid curve (e.g., by finding a local minimum or local maximum of tlie second derivative of a 
representation of the curve, such as the four*parameter logistic function). At 320, the other bound 
(e.g., endpoint) of the range is found via the second derivative of the sigmoid curve (e.g., by finding a 
10 local maximum or local maximum, inflection of ibe second derivative of a representation of tlie 
curve, such as the four^parameter logistic function). 

The usable portion, then, can be determined as the portion of the sigmoid curve between the 
two bounds (e.g., endpoints). 



15 Example 3 - Exemplary Implementation of Determination of Usable Portion 

FIG. 4 shows an excnplary inq>Iementation of determining a usable portion of a sigmoid 
curve, such as the method 300 shown in FIG. 3. In the exan^le, a standard sigmoid curve 410 has 
been generated by finding appropriate parameters for Fcxrmula (1). Such an approach is sometimes 
called "fitting a curve to the data." A second derivative 420 of the curve 410 is also shown. 

20 A first point 430 on tbe second derivative 420 mdicates a point on the curve 410 designated 

as a bound (e.g., a minimum of the second derivative 420 in the example). A second point 440 on tbe 
second derivative 430 indicates another point on the curve 410 designated as a bound (e,g., a 
maximum of the second derivative 420 in the exan^le). 

The range 460 of the curve 410 between the two bounds 430, 440 is sometimes called the 

25 "quantification range** because it is the range of the curve that is considered usable v/htn determining 
Ae concentration of a substance in a test sample. 

The point value 450 associated mfk &e point 440 is sometimes called the "quantification 
titer" because it indicates a value beyond which the curve 440 is no longer useful (e.g., a lower bound 
of concentration reliably calculable via die curve). 

30 

Example 4 - Exemplary Determination of a Threshold Value for Sigmoid Curve 
A threshold value useful in calculating concentration of a substance in a test sample can be 
detemuned in a variety of ways. One way involves finding a threshold value at which a first 
derivative for the curve reaches a benchmaric value. Such a benchmark value can be enqnrically 
35 deteimmed (e.g., chosen by a researcher based on evaluation of a number of standaid sigmoid curves) 
and then used uniformly throughout Jn one inq>lementation» the benchmark value is designated as 
&irty percent (30%) of tbe maximum of the first derivative; however, a value aj^ximatmg thirty 
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percent or a different value may be useful in other scenarios. Any of the methods can be 
iraplemented in software. 

Hie point at ^ch the first derivative is eqiial to the threshold is sometimes caUed the 
"dircshold titer*' becatise It indicates a mminniT^ detectable concentration according to the technique. 
5 Thus> if an observation of a test sample indicates a value above the threshold titer, presence of the 
substance is indicated, even if die exact concentration may not be reliably determined. Under some 
scenarios (e.g., detection of toxins), such a threshold titer can be used to advantage. 

FIG. 5 shov/s an exemplary iixiplementation of a method for finding a threshold titer 550 via 
a first derivative 520 of the related standard sigmoid curve 510. The exan:5>le portrays a scenario 
1 0 involving measurements of optical density ("OD") for various dilutions of a test s^xapk, but can be 
applied to other scenarios involving measuring concentration for a test san^>le. 

In the example, a benchmark value 540 is used to detennine a point 530 on the first 
derivative 520 at which the standard curve 5 10 crosses the benchmark value 540. From the point 530, 
a threshold titer value 550 can be detennined. If a sample has an observation above the threshold 
15 titer value 550, presence of the substance belag measin:ed is indicated 

In the example, the first derivative 520 has two points at M^ch the curve 510 crosses the 
benchmark value 540. The point related to fbs lowest concentration (e.g., the lower value of the 
curve 510) is used- 

The threshold titer value 550 can be used in a method to determine whether the substance of 
20 interest is present For exaxople, a sample for which at least one observation indicates at least such a 
vahie can be designated as containing the substance of interest Although the exanq>le describes a 
dneshold titer, the technology can be applied to metrics other &an titers. 

If desired, another threshold at 560 can be found via the first derivative and used in certaxa 
scenarios. 

25 

Example 5 - Exemplary Method Using Derivatives to Determine Concentration 

HG. 6 shows an exein|)laiy oieOu)d 600 usmg derivatives to determine subst 
concentration for a test sanqile. The me&od 600 can be in9)lementBd in software. 

At 610, a usable portion of a sigmoid curve is determined (e.g., via any of the technologies 
30 described herein). 

At 620, a tiureshold value (e.g., threshold titer) is determined (e.g., via any of the 
technologies described herein). 

At 630, for one or more observations of a saiople having unknown substance concentration, 
the observations are con9>ared to the threshold and a concentration is calculated for them via tiie 
35 usable portion of the curve. Those observations associated with a non-usable portion of the curve are 
not included. However, an observation above tire threshold can indicate presence of the substance, 
even if a concentration cannot be calculated for the observation. 
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In practice, 610 and 620 niny be pcrfonned once for a standard curve, and 630 may be 
performed plural times for plural samples (e.g., al a different tinw, location, or by a different 
computer user ai a different computer). Alternatively, in other scenarios, 610 and 620 may be 
performed plural times for the standard curve (e.g.. if the standard test serum is included plural times 
5 on plural respective plates). 

Eisample 6 - Exemplary Illustration of Back-Calculation via Sigmoid Curve 
FIG, 7 illustrates a back-calculation of concentration (e.g., use of the inverse of the four- 
pammeter logistic function) via a standard sigmoid curve. In the example, concentration of anthrax 
1 0 anti-Protective Antigen (anti-P A) IgG is being calculated to determine the concentration or presence 
of anthrax. However, the technique can be applied to any number of other scenarios for calculating 
concenti'ation of a substance. 

In the exan^le, a usable portion of the curve has been determined. For observations along 
the test-sample curve, concentrations are found. Certain observations (e.g., OD7) have values not 
15 within the usable portion of tlic curve; therefore they are not used (e.g., are discarded). Thus, a subset 
of the observations (e.g., those associated with the usable portion of the standard curve) is used to 
calculate concentratioiL 

Example 7 - Exemplary Modding of Sigmoid Curve in Optical Density Scenarios 
20 In any of the exanq)les described herem, die sigmoid curve can represent observed optical 

density for a plurality of serial dihitions. For exanqple, a reference saix^le having a known 
concentration of a substance can be observed by measuring the optical density of a plurality of 
observations for a plurality of dilutions of the sample. The sanqple can be subjected to various 
conditions before or after dilution. 
25 In some scenarios, the dilution lessens the amount of protective (e.g., cell-saving or 

antibody) agent, resulting in increased cell death. Optical density measurements at an appropriate 
wavelength can be used to measure the number of living cells. 

Example 8 - Exemplary Applications in Antibody Scenarios 
30 In any of the examples described herein, the technologies can be applied to measure the 

concentration of antibody in a san^le. For exan^)le, the sample can be diluted and added to a set of 

cells in the presence of a cell toxin. The presence of antibody will protect the cells ficom the toxm; 

the concentration of living cells thus indicates concentration of die antibody. Alternative techniques, 

such as measuring binding (e.g.> via a label or odier noarker) can be used. 
35 Presence of antibody can indicate presence of a toxin. For exanple, hi one in^lementation, 

die con^ntration of anthrax is Indirectly determined by naeasuring the presence of anti-PA IgG. 

When constmcting a standard curve for such a scenario, dilution can reduce the amount of protective 

serum (e.g., a san^le containing anti-PA IgG) that tends to preserve the life of ceUs. A low optical 
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density thus mdicates marc dead cells; a high density indicates more living cells. Less protective 
serum tends to result in more dead cells and thus a lower optical density. 

A sample having low optical density below the threshold titer is considered to test negative 
for anthrax (c,g., there was no anti-PA IgG to protect &e cells). Various other procedures can be 
5 combined into a standard protocol to promote uniformity of results across test samples. 

Esample 9 - Exemplary Technique for Determining Concentration of Antibody 
Any of the technologies described herein can be used to test for the concentration or 
presence of antibody. For cxanq>le, a method for testing a blood serum sample for the concentration 
10 of antibody in the blood sample can proceed as follows: 

(1) measuring tibe optical density of a test sani^le, wherein the test 5anQ)le is generated by 
adding the serum to cells and a toxin neutralized by the antibody, 

(2) determining whether a measurement indicating concentration of live cells {e.g., the 
optical density) fells within a usable portion of a standard sigmoid curve representing observations 

1 5 taken of a sample having a known concentration of antibody; and 

(3) responsive to deterinining the measurement indicating concentration of live celb (e.g., 
optical deixsity) feUs withia tiie usable portion, calculating a concentration via the standard sigmoid 
curve. 

In any of the examples, the antibody can be anti-PA IgG, which indicates infection by the 
20 anthrax toxin. 

During the method, one or more observations havuig an optical density outside the usable 
portion of the standard sigmoid curve can be discarded 

The usable portion of the sigmoid curve can be determined via a second derivative of the sigmoid 
curve as described in any of die exaxi:9>les. 

25 

Example 10 - Exemplary Applications 
The technologies described herein can be used m conjunction with any dilution-range-type 
micro-biological assay, including but not limited to enzyme-linked immunosoibent assay ("ELISA**), 
toxin neutralization/neutralizing assay ("TNA"), and flow opsono techniques. 
30 The technologies described herein can be applied in a variety of other bioassay analysis 

scenarios to advantage. For exan^le, the technologies can also be applied with advantage to 
polymerase chain reactive ("PGR") assay, radioimmuno assay CTOA"), cell growth assay, cell death 
assay, and flow cytometry scenarios, or any test in which sample results are derived from a standard 
sigmoid curve. 

35 
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Example 11 - Ejeniplary Ob&ervatiotis 
Observations can be taken in a variety of ways to determine concentration of a substance. 
For example, in biological scenarios, Uie number of living cells can be measured (e.g., by measuring 
optica] density at on appropriate v^avelength). 
5 AJternativeiy, a binding agent can be used. Such a binding agent may include a marker 

(e.g., chemical, fluorescent, radioactive, and the like). The marker can assist in detection. After 
binding, an observation can be taken by measuring the resulting bound substance (e.g,, by measuring 
optica] deasity). 

1 0 Example 12- Exemplary Results 

Semm samples were analyzed during 2001, 2002, and 2003 according to a protocol for a 
TNA for anti-PA antibodies. The protocol is described in Quinn et al., "Specific, sensitive, and 
quantitative enzyme-linked immunosorbent assay for human immunoglobulin G antibodies to anthrax 
toxin protective antigen". Emerging Infectious Diseases, 8:10, Centers for Disease Control, Oct 2002, 

1 5 (available at httpy/www.cdc.gov/ncidod/EID/voi8iio 1 0/02-03 80.htm) which is hereby incorporated 
herein by reference. A standard human serum was tested for neutralization of PA toxin separately on 
103 plates. 

Four samples including a serum standard (AVR414) sample were analyzed on each plate. 
Each sample consisted of a seven-point dilution with td^licates at each dilution. For samples which 

20 met primary quality control ('*QC^» four-parameter logistic curves {see Formula 1) were fitted to the 
raw data, using the NLIN procedmie of SAS software (SAS Institute, Gary NC, USA). FIG. 8 shows 
an exemplary curve 810 fit to data and flic associated confidence bands 820A, 820B. 
From tiie fitted curves, for a standard human serum tested for neutralization of PA toxin separately on 
103 plates, four titers were computed: effective-dose 50% ("ED50"), minimum detectable 

25 concentration ("MDC") and reliable detection limit ("RDL") fi:om O'Connell, the threshold titer 
described herein, the quantification titer described herein, and antibody concentration. 

The statistical observations and analytical conclusions are based on empirical data 
developed by Center for Disease Control's anthrax serology lab over the last two years. Numeric 
results apply specifically to subsets of ai^ximately 1500 anthrax TNA plates which have been 

30 analyzed in the lab and whose lab results have been analyzed statistically. 

Quantification Titer, and Threshold Titer involved less variability as measured by 
Coefficient of Variation ("CV**) and by range of values over the 103 plates. The results are shown in 
Table 1 . Titer indicates some measure of concentratioa For the nominal lower limit, MDC ranged 
from 1600 to 1 1,650 with an overall CV of 41%, v^mas the Threshold Titer ranged fiom 1690 to 

35 8,810 with a CV of 37%. For the reliable lower lunit, RDL ranged firom 1 180 to 7»660 wiOi a CV of 
38%, vMk tibe Quantification Titer ranged fi:om 1 130 to 3,960 with a CV of 29%. 
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Table I - Comparison of Iv^TOC and RDL to TT and QT (n 103) 





MDC 


TT 


RDL 


QT 


Range 


1,600-11,650 


1,690- 


1,180^ 


1,130- 




3.S10 


7,660 


3,960 


CV 


41% 


37% 


38% 


29% 



fizample 13 - Exemplary Explanation 
The described tedmologies can take advantage of the &ct that the knees (i.e., the i^per and 
5 lower bends) of a symmetric sigmoid curve (e.g.» the four-parameter logistic-log curve) can be 
identified directly through Ihe use of the first and/or second derivatives. Such an approach can be 
superior to the O'Connell approach, which evaluates error in the inverse fimction. 

The described Quantification Titer and Threshold Titer are especially precise for assays 
whose curve-fit is especially good (e.g., curves whose R-squarcs [R-squared for the non-linear curve 
10 fit is defined as the difference between the sum of squares for the observed values and the sum of 
squares for the residuals divided by the sum of squares for tiie observed values: (SSQobservBd 
SSQnsidttat) / SSQobsmd] aie in thc range of 0.99 to 1.00). 

Example 14 - Exemplary ED50, Titers, and Concentration 

15 In general, titers are computed by identitying some response threshold. Then, the fitted 

curve is used to translate that threshold firom the response axis (y or vertical axis) to the dilution axis 
(x or horizontal axis). An example is shown in FIG. 9. 

Methods of determining that response threshold are described herein. Also, the four- 
parameter logistic ("4PL'') curve was used to 1)ack-con^>ute'* (e.g., apply the inverse 4PL fimction to 

20 compute) saniple concentrations, by applying fiie relationship between response variable (optical 

density) and dilution ftom tiie standard-serum curve to Ihe observed san^le curve or raw data points. 
An example is shown in FIG. 10. 

Example 15 > Exemplary Combination of Technologies 
25 Hie described technologies can be incorporated into a TNA system. The system can 

conqmte the 4PL curve, the £D50, two in^licit titers (e.g., threshold and quantification), the usable 
portion of &e curve or quantification range^ and the in^>licit concoitmtions of test saii:ples. The 
method system can identify tiie **knees" of the 4PL curve directly througih the use of the derivatives. 
The Threshold Titer can be defined using an en^irically derived threshold in the fixst 
30 derivative to define &e point on the curve, thereby identifying a titer point on the dilution scale (e.g., 
as shown in FIG. 5). The second-derivative characteristics of the curve finther identify the knees by 
the relatively rapid change in slope. The domain between the knees is the quantification range of die 
curve. The Quantification Titer can be the dilution corresponding to Ihe upper (lower dilution) end of 
the quantification range (e.g., as shown in FIG. 4) 



35 
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Example: 16- Exemplar}' Information Regarding Using 
Sigmoid Curves to Calculate Concentration 
The fitted curve for the standaid serum, having passed primary and secondary quality 
control CQC% can establish the «or/;jo//z«</ relationship betv^reen optical density TOD") and dilution 
5 for a specific plate. Mathematically, tiiis establishes OD as a function of dilution, i.e. OD«y{DIL), 
vvherc/is the four-parameter logistic function. This is called the standard curve or characteristic 
curve for the plate. 

Further, the usable portion of tiie standard curve is that domain of dilutions or thai 
corresponding range of OD's over which the inverse ftinction can be used reliably, where "reliably" 

1 0 is defined by a specified level of precision, typically stated as a specified percent-CV across the 

implicit concentrations either across each and every dilution ("within-dilution**) or across all dilutions 
C*wjthin assay,") That is. the inverse function, DIL»/"*(OD) can be and is applied to the OD's of the 
test samples to "back-calculate" the dilutions and therefore the initial sample concentrations, since the 
starting titration dilutions are known. 

^G. 7 shows an exemplary sigmoid curve employed in a back calculation technique. A 
possible procedure for use with ti)e curve is: (1) fit the standard curve, pass QC tests, (2) define the 
usable portion of the curve, (3) for test OD's within the usable portion of the standard curve, back- 
calculate the dilutions, and (4) adjust the back-computed dilutions for the starting titration dilutiou to 
calculate the Anti-PA IgG concentrations of the test samples. Details vary acconiing to study-specific 

20 protocols. 



Example 17 - Exemplary Protocol 
Any of the examples described herein can be used m any number of protocols. For 
example, a certain number of observations falling within the usable portion of the curve can be 
25 required, various quality assurance (*'QA") tests can be required, etc. 

Example 18 - Exemplary Metrics 
In any of the examples herein, a number of metrips can be used to detennine the 
concentration of a substance. Some of the examples use a metric of optical density, but the 
30 technologies described herein can be applied to other metrics. 

Example 19 • Exemplary Indirect Measurement 
In any of the examples herein, the technologies can measure the concentration of caie 
substance to mdicate tiic concentration of another substance. Thus, measurement of a concentration 
35 can be achieved mdirectly. For example, for an antig^ such as a toxm, the cono^tration of antibody 
present in a sample taken from a subject (e.g., a human patient) can mdicate the concentration of the 
toxin in the subject's body, for example the concentration in a body fluid, such as blood, for example, 
blood serum. Further, the concentration of live cells in the presence of a toxm can indicate the 
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presence of an antibody neutralizing the toxin. Thus, the concentration of live cells in a toxin 
neutralization assay can be used to indirectly indicate the concentration of the toxia 

Esample 20 - Exemplary Approaches Using Sample Curve 
5 For any of the techiKslogies described herein, a plurality of points can be used to generate a 

test (e,g., for a $an[]|»le) sigmoid curve (e.g., according to the 4?L fcmxmla). The sample curve can 
then be used to calculate concentration in light of a usable portion of a standard sigmoid curve (e«g., 
determined via any of the technologies described herein). Such an approach can be used instead of or 
in addition to using points to calculate concentratioiL 
1 0 Although examples herein describe arrangements in \iiuch the usable portion of the standard 

curve is detemuned, an approach could alternatively or in addition determine Has usable portion of a 
test curve. 

When determining concentration in light of ^ two curves, calculations can proceed by 
considering a number of points on the test curve within tibe usable portion of the standard carve (e.g., 
1 5 using a particular point resolution). 

Even though a curve can be generated from test points filing outside the usable portion of 
the standard curve, it may be desirable to require one or more test points to fall within the usable 
portion of the standard curve (e.g., for quality control purposes). 

20 Example 21 - Exenq^lary Advantages 

Various advantages can be gained by using the exan^yles described herexa The exan4>les 
described herein may benefit from zero or more of the following advantages. 

When the described technology is used across plates, it establishes very specific quantitative 
characteristics of a given material such as a serum standard. 
25 Hie described technology allows great flexibility in determining if and how a reagent fits the 

expected model and therefore provides more robust conqnitation of endpoints. 

The described technology provides a strong quality control test for known material or 
substance. 

The described technology increases accuracy. 
30 The described technology increases repeatability. 

Given an assay which is sufficiently precise biologically, the described tedinology presented 
improves on the reproducibility and dierefore flae accuracy and precision of prior conq>utational 
methods. 

For example, v/hen the described technologies relating to first and second derivatives are 
35 applied to an anthrax TNA scenario, the long-run e£fect is that the stability of the described 
technology is higher (e.g., as measured by higher variance) tiian O'Connell's MDC and RDL. 

The described technology provides a contuiuous, as opposed to discreet en(^oint 
determination and is therefore more precise. For purposes of illustration, assume a covapzttd %CV 
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technique can only use only specific, observed dihitiom of the test-sample for end-poiutu (e.g., of tlic 
usable poition of the curve). In a worsl-casc exan^lc, if the true (but yet-unknown) usable portion of 
a curve covered the domain from a dilution of for example, 500 to 1500, and the dilution sequence for 
the test sample were 200, 400. 800. 1600.., tlicn the compared %CV approach would fail at 400 and 
fail at 1 600, leaving only the test triplicate at 800 for the back-calculadon. By many protocols, use of 
a single triplicate would be insufficient, so in this specific example, there would he no reportable test 
result. 

By contrast, tlie described technology can avoid the above outcome if the test-sample itself 
produced an acceptable sigmoid curve, in such a case, additional OD points on the test curve 
between dilutions of 500 and 799 and between 801 and 1500 can be in^utcd, using knowledge of the 
test curve gained from the dilutions 400, 1600 and beyond, even though these data are outside the 
usable portion of the standard curve. These additional points can be used to back-calculate additional 
concentrations and produce reportable results. These results rely in a predictable way on information 
gained from the test curve and from the standard curve without constraint to the two or three test 
15 dilutions points (one point in the above exanq^le) which might happen to fall in the usable domain of 
the standard curve. 

The described technology can define the usable portion of the curve before the application of 
fee curve to the test-sanple results, rather tban after the {plication of the inverse ftmction. In a 
comparison technique, 4ie measure of the usable portion of fee standard curve is whefeer or not it did 
work as evidenced by post-use error, (i.e. %CV in fee inqilicit test-sample concentrations). By 
contrast, the described technology can directly measures fee slope of fee curve and fee second- 
derivative characteristics of that slope, fecreby allowing dkect identification of fee knees of fee curve 
and fee domain between fee knees. Thus» fee usable portion of fee curve can be defined. The 
described technology can e^qploit fee characteristics of the standard curve itself, rafeer than implicit 
25 characteristics based on any particular test sample. 

Example 22 - Exemplary Concentration Units 
Alfeough fee described technologies can be used wife any concentration imits, exemplary 
units used in association wife determining concentration include titers, grams/liter (e,g., 
30 micrograms/milliiiter), or ofeer measurements of concentration. 

Example 23 - Exemplary Implementations in Software 
Any of fee mefeods described herein can be in^lemented in software. For any of fee 
methods in^lemented in software, the depicted actions can be achieved via con^uter-executable 
35 instructions carrying out fee depicted actions. Such instructions can be encoded on one or more 
conq)uter-readable media (e.g., RAM, ROM, CD-ROM, DVD, hard disk, removable media, flash 
media, and fee like). 



20 



wo 2004/084708 



FCT/US2004/0085<i6 



13 

When implemented in software, the various inputs to the methods (e.g., observations of a 
test san^le) can be provided via a user interface or sent directly via test equipment to the computer 
system. The outputs can be provided via user interface^ sent directly to output, or communicated via 
a network to another coirq>uter system, as desired. 

5 

Example 24 - Exemplary System 
FIG. 1 1 shows a system 1 100 for cairying out the technologies described herein. In the 

example, a sigmoid curve is represented by a representation 1110 (e.g., a data stmcture or set of 

variables or fields). The usable portion of the sigmoid curve is represented by a representation 1 120 
10 (e.g., a data stmcture or set of variables or fields). One or more observations of a test sample can be 

received and stored in the storage 1 130 (e.g., a database, array, or other data structure), 

A comparer 1 140 is operable to determine i^diefher an observation of 1 130 is within the 

usable portion of the sigmoid curve as indicated by the representation 1120. For exan^jle, a function 

or other logic can make such a detemiination. 
15 A calculator 1150 is operable to calculate a concentration for the test sample via the 

observations 1 1 30 and the sigmoid curve representation 1110. The calculator 1 1 50 can be configured 

so that calcvilation is performed responsive to determining that the observation is within the usable 

portion of the sigmoid cm^e. Observations determined to be outside the usable portion of the 

sigmoid curve can be rejected 
20 The usable portion of fee sigmoid carve as represented by die representation 1 120 can be 

calculated via a second derivative of the sigmoid curve as described in any of tl^ exanqiks herein. 

The system 1 100 can further include means for calcoIatEag such usable portion (e.g., by determining 

en(^ints). 

The pictured system 1 100 can be in^plemented in software stored on one or more coix:q;}uter* 
25 readable media and can be implemented in any number of computer languages. 

Example 25 - Exemplary User Interface 
FIG. 12 shows a screen shot of an exemplary user interfiice 1200 by wMch results of the 
technologies described herein can be presented (e.g., in a window 1230). 
30 In the example, analysis has been performed for six (6) samples. The window 1230 

indicates the sanople identifier, a conc^tration of a substance widiin the san:q>le, an indication of 
whether the substance was present in the 8an:q)le, the number of observations for the sample, and the 
nuniber of obsetvations discarded (e.g., for felling outside of the usable portion of die standard 
sigmoid curve). 

35 Additional, fewer, or different elements can be added to the user interfece 1200. For 

exanq>le, a representation (e.g., file name or list of parameters) of the sigmoid curve can be presented. 
Also, other discarded obserrations (e.g., filing to meet can be noted. The pr^ence of the 
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substance can be alternatively indicated (e.g.. by color, flashing text, or icon). Similarly, the otlicr 
results can be alternatively indicated. 

Ezaniple 26 - Other Esemplary User Interfaces 
5 FIG, 1 3 shows a screen shot of an exemplary user interface 1 300 for presenting results 

related to a standard curve. In tlie example, the usable portion of the curve between 1310 and 1320 is 
depicted in one color (e.g., blue), and the points outside the useable portion but within the thresholds 
are depicted in another color (e.g., brown). The remainder of tlie curve is in yet another color (e.g., 
black). 

10 FIG. 14 shows a screen shot of an exemplary user interface 1400 similar to that of FIG. 13, 

but for a test curve. Again, the usable portion of the curve betweejD 1410 and 1420 can be depicted in 
one color (e.g.» blueX and the points outside Ute usable portion but within Hbc thresholds can be 
depicted in another color (e.g., brown). The remainder of tbe curve can be in yet another color (e.g., 
black). 

15 FIGS. 15 and 16 show a screen shot of an exemplary user interface 1 500, 1 600 that presents 

results in text form. The interface includes both point-based IgG concentration (e.g., determined 
applying actual observations to a standard sigmoid curve) and curve-based IgG concentration (e.g., 
determined applying a curve fit to actual observations to a standard sigmoid curve). 

Any combination of the user inter&ce described herein can be used in software. Additional, 

20 fewer, or dififerent elements can be used as an alternative to those pictured 

Example 27 - Exemplary Alternatives 

For those actions specified as coirputer'executable, such actions can be performed fully- 
automatically (e.g., without human intervention) or semi-automatically (e.g., with assistance from a 
25 human operator). One or more con^uter-readable media can comprise the instructions described as 
con^uter-executable. 

Various in^lementations of tiie technologies described herem can be called the ''Taylor 
method,** named after inventor Thomas H. Taylor, Jr. 

in view of the many possible embodhxients to which the principles of the invention may be 
30 appUed, it should be recognized that the Olustrated embodinoents are examples of the invention, and 
should not be taken as a limitation on the scope of the invention. Rather, the scope of the invention 
includes what is covered by the following clainos. I therefore claim as my invention all that comes 
within the scope and spirit of these claims. 
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We claim: 

1. One or more computer-readable media comprising conputer-executable 
instructions for pcrfonning a method to calculate concentration of a substance in a test sample, the 
method comprising: 

for at least one observation of a metric for ^ test sample, finding where on a usable portion 
of a standard sigmoid curve the observation lies, wherein ^e usable portion of the standard sigmoid 
curve is determined via a second derivative of the standard sigmoid curve; and 

based on a location of the observation on the standard sigmoid curve, calculating a 
concentration of the substance. 

2. The one or more con^iuter-readable media of claim 1 wherein it& sigmK>id curve is 
represented via a four-parameter formula. 

3. 'fhe one or more con^uter-readable media of claim 1 whereia the standard sigmoid 
15 curve represents a sigmoid curve fit to a plurality of observations taken of a refi^ence sample having 

a known concentration of the substance. 

4. The one or more conq}ater-readable media of claim 1 further comprising con^)uter- 
executable instructions forperfonning the following: 

20 determining for at least one observation of a metric for the test san^ile whether tiie 

observation is above a threshold vahie^ wherein the tixreshold value is determined via a first derivative 
of the standard sigmoid curve; and 

indicating whether the obs^ation is above the threshold vahie. 

25 5. The one or more con:q)uter-'readable media of claim 1 wherein: 

the observation indicates optical density for the test sample. 

6. The one or more computer-readable media of claim 5 wherein: 
the concentration indicates an amount of antibody in the test sanqile. 

30 

7. The one or more con^uter-readable n^dia of claim 6 wherein: 
the concentration radicates an amount of anti-PA IgG in the test sample. 

8. One or more con^uter-readable media comprising conputer-executable 

35 instructions for performing a method to calculate concentration of a substance in a test sairqple, the 
method comprising: 

for a plurality of observations of a metric for the test sanQ>le« fitting a test sigmoid curve to 
the observations; and 
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calculating a concentration of tlic substance in the test sample via tlie lest sigmoid curve and 
a usable portion of a standard curve, wherein the usable portion of the standard sigmoid curve is 
detemtined via a second derivative of the standard sigmoid curve. 

5 9. The one or more con^utcr-rcadable media of claim 8 ftirther comtphsing computer- 

executable instructions for performing the following: 
indicating tlie concentration of the substance. 

1 0. The one or more computer-readable media of claim 8 further comprising computer- 
3 0 executable instructions for performing the following: 

displaying the concentration of the substance. 

1 1 . One or more computer-readable media con^rising connputer-executable 
instructions for performing a method to calculate concentration of a substance in a test sample, the 
method con^rising: 

finding a usable portion of a sigmoid curve, wherein the usable portion of the sigmoid curve 
is determined via a second derivative of the sigmoid curve; and 

calculating a concentration of the substance in the test sanaple via the usable portion of the 
sigmoid curve, 

12. One or more conqjuter-ieadable media con^rising computer-executable 
instructions for performing a method con^nsing: 

for a plurality of dilutions of a test sample, receiving respective measurements of optical 
density indicating concentration of live cells within the dilutions; 

via the measurements, calculating a concentration of anti-PA IgG for die test san^le via a 
usable portion of a sigmoid curve representing concentrations of live cells within dilutions of a 
reference sample havmg a known quantity of anti-PA IgG, wherein the sigmoid curve is represented 
via a four-parameter logistic technique, and wherein a usable portion of die sigmoid curve is 
determined via a second derivative of the sigmoid curve; and 

indicating the concentration of anti-?A IgG for the test sample. 

13. A computer-implemented method of calculating concentration of a substance in a 
test sanple having an unknown concentration of the substance, flie method comprising: 

detennining a usable portion of a sigmoid curve fit to data points representing observations 
35 of a refer^ce san^le having a known ccmcentration of the substance; and 

calculating the concentration of the substance in fhs test sample based on a subset of 
observations of the test sample, wherein die subset is associated witib die usable portion of die 
sigmoid curve. 



20 



25 
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14. The method of claim 13 further con^jrising; 

excluding at least one excluded observation of ti^ic test sample responsive to detenninijog the 
excluded observation is outside the usable portion of the sigmoid curve. 

5 

15. The method of claim 13 wherein determining a usable portion of the sigmoid curve 
comprises calculating a second derivative for the sigmoid curve. 

16. The method of claim 13 wherein determining a usable portion of the sigmoid curve 
1 0 corr^)rises designating a portion between a minimum and a maximum of a second derivative for the 

sigmoid curve as the usable portion of the sigmoid curve. 

17. The method of claim 13 wherein a point on the*8igmoid curve relatrng to a 
threshold for a first derivative of the sigmoid curve is used as a lower threshold to indicate presence 

15 of the substance. 

18. A con^^uter-Lo^lemented mediod of determining the concentration of antibody in a 
blood serum sanq)Ie, the method con^msxng: 

receiving a measurement of concentration of live cells in a test sample, wherein the test 
20 sample is ger^ted by adding the serum to cells and a toxin neutralized by the ant&od;^ 

detennining whe^r the concentration of live cells &lls within a usable portion of a standard 
sigmoid curve representing observations taken of a sample having a known concentration of 
antibody; and 

responsive to determining the concentration of live cells falls within the usable portion, 
25 calculating a concentration via the standard sigmoid curve. 

19. One or more conqmter-readable media having C0n9)uter-executable instructions for 
performing the method of claim 18. 

30 20. The method of claim 1 8 wherein results for plural test samples for phiral dilutions 

of an original test sample are included in the calculating. 

21. The method of claim 18 vi^erein concentration of live cells is indicated by optical 
density of die test san:q>le. 

35 

22. The method of claim 18 wherein concentration of live cells is indicated by optical 
density of the test sanq>le. 
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23. The method of cJatm 1 8 wherein the antibody is anti-P A IgG. 

24. The method of claim 1 8 further comprising: 

discarding at least one observatloQ having a concentration of live cells outside tlie usable 
5 portion of the standard sigmoid curve. 

25. Tlie method of claim 1 8 further comprising: 

in software, determining the usable portion of the sigmoid curve via a second derivative of 
the sigmoid curve. 

10 

26. A sofhvare system encoded on one or more con^uter-readable media, the software 
system comprising: 

a representation of a characteristic sigmoid curve; 

means for designating the usable portion of the characteristic sigmoid curve; 
1 5 means for receiving at least one observation of a test sample; 

means for determining whether the observation of the test sample is within the usable 
portion of the duiracteristic sigmoid curve; and 

means for calculating a concentration for the observation responsive to determining that the 
observation is within the usable portion of dse characteristic sigmoid curve. 

20 

27. The software system of claim 26 wherein the usable portion of the characteristic 
curve is calculated via a second derivative of the sigmoid curve. 

28. The software system of claim 26 further coniiprising: 

2^ means for determining the usable portion of the sigmoid curve via a second derivative of die 

sigmoid curve. 



29. The software system of claim 26 further comprising: 

txseans for rejecting an observation responsive to determining ^t the observation is outside 
30 the usable portion of the characteristic sigmoid curve. 
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30. One or more conputer-readabJe media con^nising computer-executable 
instructions for performing a method to indicate presence of a substance in a test san^ie, the method 
comprising: 

for at lea3t one observation of a metric for the test sample, detennining wheflier the 
5 observation is higher than a threshold value, wherein the threshold value is detennined via a first 
derivative of a standard sigmoid curve; and 

responsive to determining the observation is higher than the threshold value, indicating 
presence of the substance. 
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